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Abstract. We use a rigorous method, based on Bayesian infer- 
ence, for calculating the odds favoring the hypothesis that any 
particular class of astronomical transients produce gamma-ray 
bursts over the hypothesis that they do not. We then apply this 
method to a sample of 83 Type la supernovae and a sample of 
20 Type Ib-Ic supernovae. We find overwhelming odds against 
the hypothesis that all Type la supernovae produce gamma- 
ray bursts, whether at low redshift (10^ : 1) or high-redshift 
(10^^ : 1), and very large odds (6000 : 1) against the hypoth- 
esis that all Type lb, Ib/c, and Ic supernovae produce observ- 
able gamma-ray bursts. We find large odds (34 : 1) against the 
hypothesis that a fraction of Type la supernovae produce ob- 
servable gamma-ray bursts, and moderate odds (6 : 1) against 
the hypothesis that a fraction of Type Ib-Ic supernovae produce 
observable bursts. 
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1. Introduction 



The discovery of supernova 1998bw (Galama et al. 1998) 
within the 8 arcminute radius of the BepppoSAX WFC error 
circle for GRB 980425 (Soffitta et al. |1998[ ) has led to the hy- 
pothesis that some GRB sources are Type Ib-Ic SNe. 

There are some serious difficulties with this interpretation 
of the data for GRB 980425/SN 1998bw: the supernova oc- 
curred outside the NFI error circle of a fading X-ray source 



(Pian et al. 1998a 1998b 1999 ; Piro et al.l998). This source 



had a temporal decay consistent with a power-law index of 



1.2 (Pian et al. 1998b), which resembles the temporal be- 



havior of X-ray afterglows seen in almost every other GRB fol- 
lowed up with the SAX NFI. It must therefore be viewed as a 
strong candidate to be the X-ray afterglow of GRB 980425. 

Moreover, if the association between GRB 980425 and SN 
1998bw were true, the luminosity of this burst would be ~ lO''^^ 
erg s~^ and its energy would be ^ 10**^ erg. Each would there- 
fore be five orders of magnitude less than that of other bursts, 
and the behavior of the X-ray and optical afterglow would be 
very different from those of the other BeppoSAX bursts, yet 
the burst itself is indistinguishable from other BeppoSAX and 



BATSE GRBs with respect to duration, time history, spectral 
shape, peak flux, and and fluence (Galama et al. 1998). 

In view of these difficulties, the safest procedure is to regard 
the association as a hypothesis that is to be tested by search- 
ing for correlations between SNe and GRB in catalogs of SNe 
and GRB s, exclu ding SN 1998bw and GR B 9804 25. Wang & 
Wheeler (1998a; see also Wang & Wheeler 1998b) have corre- 
lated BATSE GRB with Type la and with Type Ib-Ic Sne. They 
found that the data was "consistent" with the assumption of an 
association between GRB and Type Ib-Ic SNe. 

In the present work we improve upon the Wang & Wheeler 
correlative study by introducing an analysis method based on 
Bayesian inference, and therefore using the likelihood function, 
that incorporates information about the BATSE position errors 
in a non-arbitrary way and that is free of the ambiguities of a 
posteriori statistics. The method also accounts the fact that the 
BATSE temporal exposure is less than unity. 

2. Methodology 



We use the BATSE 4B GRB Catalog (Meegan et al. |1998[ ), 
and BATSE bursts that occurred subsequent to the 4B cata- 
log but before 1 May 1998. We also use the Ulysses supple- 
ment to the BATSE 4B catalog, which contains 219 BATSE 
burst s for w hich 3rd IPN annuli have been determined (Hurley 
et al. 1998[ ). Hurley (private communication, 1998) has kindly 
made available at our request 3rd IPN annuli for an additional 
9 BATSE bursts that occurred subsequent to the period of the 
BATSE 4B catalog but before 1 May 1998. 

We have compiled three Type I SNe samples. The first is a 
sample of 37 Type la SNe at low redshift (z < 0.1) from the 
CfA SN Search Team (Riess 1998, private communication). 
The second is a sample of 46 moderate redshift (0.1 < z < 
0.830) Type la SNe from the Supernova Cosmology Project 
(Perlmutter 1998, private communication). The third sample 
consists of 20 Type lb, Ib/c, and Ic SNe compiled from lAU 
Circulars and various SNe catalogs. 

We compare three hypotheses: 
Hi: The association between SNe and GRBs is real. If a SN is 
observed, there is a chance e that BATSE sees the associated 
GRB, where e is the average BATSE temporal exposure. While 
e varies with Declination, the variation is modest and we ne- 
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gleet it. The probability density for the time of occurrence of 

the ith supernova is assumed uniform in the interval [ti , + r^] , 

so that all GRBs that occur in that interval have an equal prior 

probability of being associated with the SN. 

H[: The association is real, but only a fraction / of detectable 

SNe produce detectable GRB. 

H2'- There is no association between SNe and GRBs. 

We calculate the Bayesian Odds, O favoring Hi over H2- 
The details of the calculation are presented in Graziani et al 
(1998). We separately compare H[ to H2, denoting the odds 
favoring H'l over H2 by C. Finally, we calculate the posterior 
probability distribution for the parameter / in the model H[, 
and infer an upper limit on its value. 

3. Results 

We find overwhelming odds against the hypothesis that all Type 
la supernovae produce gamma-ray bursts, whether at low red- 
shift (O = 10^ : 1) or high-redshift (O = lO^^ : 1). We also 
find large odds (C = 34 : 1) against the hypothesis that a 
fraction of Type la supernovae produce observable gamma-ray 
bursts 

We find very large odds (O ~ 6000 : 1) against the hypoth- 
esis that all Type lb, Ib/c, and Ic supernovae produce observable 
gamma-ray bursts. We also find moderate odds (C = 6:1) 
against the hypothesis that a fraction of Type Ib-Ic supernovae 
produce observable bursts. If we nevertheless assume that this 
hypothesis is correct, we find that the fraction /sn of Type lb, 
Ib/c and Ic SNe that produce observable GRBs must be less 
than 0.17, 0.42, and 0.70 with 68%, 95%, and 99.6% proba- 
bility, respectively. These limits are relatively weak because of 
the modest size (20 events) of our sample of Type Ib-Ic SNe. 



The parameter / represents the fraction of observable SNe 
which can produce observable GRBs. Some observable SNe 
might not produce observable GRBs because of intrinsic ef- 
fects, such as beaming, whereas others might not do so because 
the sampling distance for SN-produced GRB could be less than 
the sampling distance for the SNe themselves. It is possible to 
separate these two effects, by writing / = /intrinsic x /sampling- 
Obviously, the more interesting quantity is /intrinsic, since it 
addresses the nature of the GRB sources. 

We may attempt to find a constraint on /intrinsic by assum- 
ing the correctness of the identification of SN 1998bw with 
GRB 980425 and using that association to estimate /sampling, 
under the (dubious) assumption that the GRBs produced by 
Type Ib-Ic SNe are standard candles. The result is /sampling ~ 
1.9 X 10^"^. This is rather bad news for the prospect of con- 
straining /intrinsic, sinCC the product of /intrinsic and /sampling 

is only constrained by the data to be less than about 0.7. Thus 
this argument can place no constraint on /intrinsic- Elsewhere, 
we show that placing a meaningful constraint on /intrinsic 
would require a GRB experiment approximately 80 times more 
sensitive than BATSE (Graziani et al. |1998|). 



One can also approach the proposed association between 
SNe and GRBs from the opposite direction. The interesting 
question, from this point of view, is what fraction /grb of 
the GRBs detected by BATSE could have been produced by 
Type Ib-Ic SNe? (note that this is different from the question 



addressed by Kippen et al. [1998], who constrained the frac- 



tion of BATSE GRB that could have been produced by known 
SNe). Such a limit may be derived under the same assumptions 
as made in the previous paragraph. We find that no more than 
90 SNe could have produced GRB detectable by BATSE, 
indicating that /grb can be no more than about 5%. 



4. Discussion 

We find large odds against the hypothesis that all Type lb, Ib/c, 
and Ic supernovae produce observable gamma-ray bur sts, the 
specific hypothesis considered by Wang & Wheeler (1998c, 
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